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Abstract. The superconducting fluctuations well inside the normal state of Fe-based 
superconductors were studied through measurements of the in-plane paraconductivity 
and magnetoconductivity in high quality BaFe2- :r Ni : jAs2 crystals with doping levels 
from the optimal (x = 0.10) up to the highly overdoped (x = 0.20). These 
measurements, performed in magnetic fields up to 9 T perpendicular to the ab (Fe) 
layers, allowed a reliable check of the applicability to iron-based superconductors of 
Ginzburg-Landau approaches for 3D anisotropic compounds, even at high reduced 
temperatures and magnetic fields. Our results also allowed us to gain valuable insights 
into the dependence on the doping level of some central superconducting parameters 
(coherence lengths and anisotropy factor). 
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1. Introduction 

As the pairing mechanism of Fe-based superconductors is not yet established, the 
phenomenological descriptions of their superconducting transition are at the forefront 
of the research in these materials. [T] A powerful tool to probe these descriptions is the 
superconducting fluctuation effects that appear in the normal state. (2j El HI El El El 
13 HOI El E2] At present, some experimental aspects of these fluctuations have been 
studied through observables like the magnetization [31 HI El E], the specific heat [7] or 
the electrical conductivity (3 El [101 HH [12]. However, their behavior at high reduced 
temperatures and magnetic fields (in the short wavelength fluctuation regime) and 
their onset temperature, remains at present unexplored at a quantitative level. The 
importance of these aspects, which include the possible presence of phase fluctuations 
well above the superconducting critical temperature (T c ),j3j HI El E] is enhanced by the 
comparison with the high-T c cuprate superconductors (HTSC), for which the onset and 
the influence of doping on their superconducting fluctuations is at present one of the 
most central and debated issues of their phenomenology. [131 EH H3 EE E] 

The central aim of this paper is to present detailed experimental results on the 
superconducting fluctuations above T c in iron-based superconductors as a function of 
the doping level. To achieve this, we have measured in an extended temperature region 
above the superconducting transition the in-plane fluctuation electric conductivity, 
Aa ab (T, H), in optimally-doped and overdoped BaFe 2 - a ;Ni a ;As 2 crystals (0.1 < x < 0.2), 
under magnetic fields up to 9 T perpendicular to the ab (Fe) layers. We have 
also used these experimental data to probe the phenomenological descriptions of the 
superconducting fluctuations based on the Gaussian Ginzburg-Landau (GGL) approach, 
adapted to the 3D anisotropic nature of these compounds, and also to take into account 
the short wavelength fluctuation regime. [HI [19] Our analysis allows us to gain valuable 
insight into the doping dependence of their superconducting parameters. In particular, 
the anisotropy factor is shown to increase significantly in strongly overdoped samples. 

2. Experimental details and results 

2.1. Crystals fabrication and characterization 

The BaFe2-xNi x As2 samples used in this work are plate-like single crystals (typically 
5 x 2 x 0.1 mm 3 ) with the c crystallographic axis perpendicular to their largest face. 
They were cleaved from larger single crystals grown by the self-flux method. Their 
nominal Ni doping levels are x = 0.10, 0.15, 0.18 and 0.20, although the real doping 
level was found to be a factor ~ 0.8 smaller (see Ref. [20], where all the details of 
the growth procedure and characterization may be found). We checked the excellent 
stoichiometric and structural quality of the crystals studied here by x-ray diffraction. In 
particular, the (00/) linewidths were found to be slightly larger [A(20) ~ 0.10°, FWHM] 
than the corresponding instrumental linewidths, A(26) ~ 0.07°. This was attributed 
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to a dispersion in the c-axis lattice parameter L C J|| and was used to roughly estimate 
Ax ~ 1CT 2 through the L c (x) dependence presented in Ref. [20]. In turn, the rocking 
curves for the (008) lines indicated that the dispersion in the c-axis orientation was 
lower than ~ 0.05°. 

2.2. Resistivity measurements: Determination of the superconducting transition 
temperatures and transition widths 

The in-plane resistivity (along the ab layers) was measured in the presence of magnetic 
fields up to 9 T perpendicular to the ab layers with a Quantum Design physical property 
measurement system (PPMS) by using four contacts with an in-line configuration and 
an excitation current of ~ 1 mA at 23 Hz. The uncertainties in the geometry and 
dimensions of the crystals, and the finite size of the electrical contacts (stripes typically 
0.5 mm wide) lead to an uncertainty in the p ab amplitude of ~ 25%. An example of 
the Pab(T) dependence around the transition (corresponding to the crystal with x=0.1) 
is presented in Fig. 1(a). The zero- field transition temperature, T c , was estimated 
from the maximum of the corresponding dp a b/dT curve, and the transition width as 
AT C m 2(T C — T~), where T~ is the temperature at which p ab = is attained. In the 
optimally doped crystal (x = 0.1) T c = 20.0 K and AT C w 0.3 K, which is among 
the best in the literature for crystals of the same composition. [221 123] EH EE] This 
allows investigation of fluctuation effects down to reduced temperatures, e = ln(T/T c ), 
of the order of AT C /T C w 10~ 2 . Overdoped crystals present a slightly wider resistive 
transition (AT C « 0.6 K) probably due to the T c (x) dependence and an x-distribution 
(taking into account that in the overdoped regime \dT c /dx\ ~ 200 K it may be estimated 
Ax ~ 0.003). Nevertheless, even in these samples, fluctuation effects may be studied in 
a wide range of temperatures above the zero-field T c by just applying a magnetic field of 
the order of ~1 T due to the T C (H) shift to lower temperatures. A detailed analysis of 
the effect on the measured fluctuation conductivity of the uncertainty in T c is presented 
in Appendix A. 

2.3. Determination of the fluctuation contribution to the electric conductivity 

An overview of the p a b{T) data up to room temperature is presented in the inset of 
Fig. 1(a). The seemingly non- monotonous x-dependence at 300 K may be explained 
in terms of the above mentioned geometrical uncertainties. As expected. |20[ [26] for 
x > 0.1, the kinks associated with structural (tetragonal-orthorhombic) and magnetic 
(paramagnetic-antiferromagnetic) transitions are not observed, and for x = 0.1 they 
are irrelevant when compared with the rounding associated to the fluctuations. This 
is an important experimental advantage in order to determine the superconducting 
contribution to the electric conductivity, Aa ab = p~ b l — p~ b B , where p a b,B is the normal- 
state or background contribution. In view of the linear temperature dependence of 

| A similar procedure was previously used to investigate compositional inhomogeneities in non- 
stoichiometric high-T c cuprates, see e.g., Ref. [21] 
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dp a b/dT (an example for x = 0.1 is presented in Fig. [1(b) ), the background contributions 
were estimated by fitting a quadratic polynomial to the measured p a b{T)u from ~ 4T C 
down to T onset , the temperature below which fluctuation effects are measurable. In turn, 
this was determined as the temperature at which dp a b/dT rises above the extrapolated 
normal-state behavior beyond the noise level [see Fig. 1(b)]. As shown in Fig. 1(c), 
T onse t{x) ~ 1.5T c (x) an issue that will be analyzed later. An example of the background 
contribution (for the x — 0.1 crystal) is shown in Fig. 1(a). On the scale of the effect 
of the superconducting fluctuations, p a b,B{T) is independent of the applied magnetic 
field. This further experimental advantage to study the fluctuation-induced in-plane 
magnetoconductivity is confirmed by the detailed representation of Fig. 2(a), where it 
is shown that a a b(0)—a a b(H) becomes negligible well above T c . A detailed analysis of the 
effect on the measured Aa a b of the typical uncertainty in the background contribution 
is presented in Appendix A. 



3. Data analysis 

3.1. Comparison with the conventional Aslamazov-Larkin approach 

An example (for the x = 0.1 crystal) of the Aa a b(H) dependence at several temperatures 
just above T c is presented in Fig. 2(b). In the H — > limit Aa a b tends to a constant 
value, in qualitative agreement with the conventional Aslamazov-Larkin (AL) result. [27] 

where e is the electron charge, fi the reduced Planck constant, and £ c (0) the c-axis 
coherence length amplitude. However, Acr ab decreases above a temperature-dependent 
magnetic field which is close to the scale for the observation of finite-field effects (see 
below) J§| An example of the Aa a b dependence on e (also corresponding to the x = 0.1 
crystal) is shown in Fig. 3(a). The double- logarithmic scale was chosen to explore in 
detail the high-e region. As may be clearly seen in this figure, the fit of Eq. (JT]) to the low- 
field data (dashed line) is excellent at low reduced temperatures (up to e ~ 0.06), except 
for e < AT C /T C « 0.015 where T c inhomogeneities may play a role (see below). However, 
a large discrepancy is found at high-e values: while the AL theory never vanishes, the 
experimental Aa a b ends up below the experimental resolution for e ~ 0.3 for all the 
applied magnetic fields. A similar rapid falloff of the fluctuation induced conductivity 
was already observed more than 30 years ago in low-T c conventional superconductors by 
Johnson and coworkers, [29] [30] who already stressed that this type of behaviour could 



§ It is worth noting that a negative fluctuation magnetoresistance has been observed in granular 
metallic films in high fields (of the order of the critical field) , see Ref. [28] . This effect is not observed here 
even in the vicinity of T c , where superconducting and normal domains may coexist as a consequence 
of T c inhomogeneities. On the one side, the magnetic field used in our experiments may not be 
strong enough as to observe such an effect. On the other, the indirect contributions to the fluctuation 
conductivity needed to explain the negative magnetoresistance (see Ref. [55]) may not be present in 
non-s wave superconductors, as it may be the case of iron pnictides. 
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not be described in terms of a power law in e. Thus, our present data provides further 
evidence of the failure of those proposals which, like the momentum cutoff approach[31J 
and equivalent microscopic calculations [32] [jj, predict a definite critical exponent for 
Aa ab at high reduced temperatures. 

3.2. Comparison with the GL theory with an energy cutoff in the fluctuations spectrum 

The failure of GL-based approaches to explain the A<r a &-data at high-e has also been 
observed in HTSC,[18] and it also appears in the fluctuation diamagnetism of both 
HTSC (33J [3l] and low-T c alloys. [351 136] It has been attributed to the fact that the 
GL-theory overestimates the contribution of the high-energy fluctuation modes. [35] In 
fact, although GL approaches are formally valid only in the vicinity of the transition, 
it was found that its applicability may be extended to the high-e region through the 
introduction of an energy cutoff. [HI [191 ESI El] Since the energy of the fluctuation 
modes increases with H, the inclusion of such a cutoff is also needed when analyzing 
the effect of a finite applied magnetic field on the superconducting fluctuations. [37] In 
the case of 3D anisotropic superconductors, the one well adapted to 122 iron arsenides, 
the paraconductivity in presence of an energy cutoff may be easily calculated by using 
the procedure proposed in the pioneering work by A. Schmid. [38] These calculations, 
whose details are presented in Appendix B, lead to 



where ip 1 is the first derivative of the digamma function, h = H/[<p / '2^0^(0)] is the 
reduced magnetic field, £ a fc(0) is the in-plane coherence length amplitude, and c is the 
cutoff constant (expected to be of the order of 0.5) [19]. In the zero magnetic field limit 
(i.e., for h <C e) and in the absence of cutoff (c — > oo), Eq. (T5]) reduces to the AL 
expression, Eq. (TT]). It is also worth noting that Eq. (TJJ) leads to the Aa ab vanishing 
at £ — c. In an attempt to check the applicability range of the 3D-anisotropic GL 
approach under an energy cutoff, in what follows we will compare our experimental 
data with Eq. 

3.2.1. Optimally doped sample: A first check may be carried out through measurements 
performed with h = 0, because in this case Eq. (T5]) depends only on £ c (0) and c. In the 

|| As we have already commented in Ref. |18j . the microscopic approach proposed in this work results 
to be equivalent to apply a momentum cutoff in to the GL-theory (see Refs. [31]), since both lead to 
an asymptotic behaviour of the paraconductivity at high reduced temperatures proportional to epsilon 
to the -3. Such a behavior was seemingly observed in cuprate superconductors by Varlamov and 
coworkers (see Fig. 7.6, and the corresponding reference, in that book). Nevertheless, such a result, 
the only accounted on the paraconductivity in that book, seems to be an artifact of the data analysis 
and, in any case, it is in contradiction with the measurements published until now by other groups in 
any low- or high-T c superconductor (for earlier results see, e.g., Ref. |31j and references therein). 




(2) 
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optimally doped sample the fit of Eq. (J2J) to the Aa ab (e, h = 0) data is excellent above 
e = 0.02, including the Acr ab vanishing at high e-values [see Fig. 3(a)]. This analysis 
leads to £ c (0) = 0.8 nm and c = 0.39. By using these values, Eq. (j2J) was then fitted to 
the data obtained with h > with £ a b(0) as the only free parameter. The fit quality is 
also excellent up to the largest field used in the experiments leading to £ a f>(0) = 2.3 nm. 
For comparison, the conventional AL approach [Eq. (JTJ)] evaluated with the same £ a b(0) 
and £ c (0) values is represented as a dotted line in Fig. 3(a). 

As expected, the adequacy of Eq. (J2J) extends to the Aa a b(H) e representation of 
Fig. 2(b), where the solid lines were evaluated by using in Eq. (j5J) the above £ a &(0), £ c (0) 
and c values. The dashed line in this figure is the crossover to the region at which finite 
field effects are expected to be relevant, and was evaluated by using h = e in Eq. It 
is worth mentioning that the in-plane magnetoconductivity, cr a b(0) — a a b(H), in Fig. 2(a) 
is also in excellent agreement with Aa a b(0) — Aa a b(H) as evaluated from Eq. (J2J) by using 
the same £ a b(0), £ c (0) and c values (solid lines). As the a a b(0) — o~ a b{H) data do not 
depend on a background subtraction, this fact further validates the procedure used to 
determine the background contribution to obtain Aa ah in Figs. 2(b) and 3. 

3.2.2. Overdoped samples: An example of the Aa a b{e)n dependence corresponding 
to an overdoped crystal (x = 0.15) is shown in Fig. 3(b). The upturns observed in 
the low- if isofields (indicated with arrows) are associated with the above mentioned 
T c inhomogeneities inherent to non-optimally-doped samples. It is expected that a 
magnetic field of the order of H c2 (0)AT c /T c > 1 T will shift T C (H) so that Aa a b would 
be unaffected by T c inhomogeneities down to e — 0. Then, we have fitted Eq. ([2]) to 
the experimental data for fi H > 3 T with £ a &(0), £ c (0), and c as free parameters. The 
agreement is excellent in the entire e range, and even extends to the lowest magnetic 
fields for e-values well above AT C /T C . A similar result is obtained in the other overdoped 
crystals studied. The x-dependence of the resulting £ a &(0), £ c (0) and c will be analyzed 
in the next Section. 

We have checked that the T c -inhomogeneities model proposed in Ref. [10], which 
is based on Bruggeman's effective medium theory, [H] accounts for the upturn in 
the H = measurement. According to this model, the temperature dependence of 
effective electrical conductivity in presence of T c inhomogeneities, a°f(T), is obtained 
by numerically solving 



where S(T C ) is the distribution of critical temperatures in the sample, and a a b{T,T c ) = 
Pab B (T) + Aa ab (T, T c ) with Acr ab (T,T c ) given by Eq. ([2]). The dot-dashed line shown in 
Fig. 3(b) was obtained by just assuming for 5(T C ) a Gaussian distribution 1.2 K wide 

% See e.g., Ref. |35) . This criterion for the presence of finite- field effects may be more intuitively written 
as H = 0o/27r/io^ab(7 1 ). This magnetic field scale is sometimes referred to as ghost critical field, see 




(3) 



e.g., Ref. [35]. 
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(FWHM), which is in reasonable agreement with the AT C value determined for this 
sample. 

4. Discussion of the results 

4-1. Dependence of the superconducting parameters on the doping level 

The x-dependence of the coherence lengths resulting from the above analysis is presented 
in the inset in Fig. 3(a). £ a b(0) grows monotonically with x, whereas £ c (0) is roughly 
independent of x. As a consequence, the anisotropy factor 7 = £ a f>(0)/£ c (0) grows 
with x from ~ 2.9 in the optimally doped crystal, up to above ~ 10 in the most 
overdoped crystal (x = 0.2). Such a large 7 value is, however, consistent with the 
3D behavior of this material because the corresponding 2£ c (0) is still larger than the 
Fe-layers separation. Previous studies in BaFe2- :r Ni r As2 focus on samples close to 
optimal doping, and report superconducting parameters consistent with our present 
results. [22l [23j EU [25] Although the large 7 value observed in the highly overdoped 
BaFe1.8Nio.2As2 was never observed in a 122 compound, some studies in the very similar 
BaFe2_a;Co x As2 reported that 7 increases on increasing x slightly above the optimal 
value for this compound. (2U 22} H3] This effect was associated with the change in the 
structural/magnetic state of the system. [24"| |4*2] or with an increase in the ratio of 
inter /intraband coupling. |43] 

4-2. Temperature onset for the fluctuation effects 

Given the relation c = ln(T onset /T c ), the cutoff constant (c ~ 0.3 — 0.7) turned 
out to be consistent with the T onset values in Fig. 1, roughly 1.5T c (x). In a recent 
study of other iron pnictides (Ref. [33]) a change is observed in the Raman modes 
below a comparable temperature (~ 2T C ) which is also attributed to superconducting 
fluctuations. Remarkably, such a T onset /T c ratio is also within the values observed in 
other superconducting families, including conventional low-T c superconductors [36] 137] . 
MgB 2 [35], NbSe2 [36], and HTSC not severely affected by T c inhomogeneities. [341 I4T] 
El This study extends to Fe-based superconductors the proposal that superconducting 
fluctuations vanish at the temperature at which the superconducting wavefuntion 
shrinks to lengths of the order of the pairs size. [19] 

4-3. Relevance of phase fluctuations or indirect contributions to Aa 

A direct consequence of the excellent agreement of the GL theory to explain our data 
is the absence of appreciable indirect contributions to the fluctuation-induced in-plane 

+ Note that some works report the observation of e onso t-values different from 0.5 in non-optimally- 
doped HTSC and amorphous low-T c superconductors. See e.g., Ref. [48]. However, these differences in 
the £ nset values could be related to the strong T c dependence on the stoichiometry in these materials, 
and the unavoidable presence of stoichiometric inhomogeneities, see Ref. [21] and [49] . 
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conductivity and magnetoconductivity of Fe-based superconductors, including the Maki- 
Thompson (MT) and the so-called density of states (DOS) contributions. These indirect 
contributions have been found to be negligible also in HTSC,[5n] and in both systems 
this could be attributed to the non s-wave pairing. [5Tj 

Finally, it has recently been claimed that conventional GL approaches for the 
fluctuation effects above T c are not applicable at low field amplitudes (typically below 
1 T) in SmFeAsOo.sFo.2, due to the presence of phase fluctuations. [6] The relevance of 
phase fluctuations has also been proposed for a member of the less anisotropic 122 family, 
Bai_ x K x Fe2As2.|3] However, our present results show that the anomalous increase of 
fluctuation effects at low field amplitudes may be explained in the framework of GL 
approaches by taking into account the effect of T c inhomogeneities. This could suggest 
that the effect of phase fluctuations in these materials is less important than previously 
estimated. 

5. Conclusions 

We have presented detailed measurements of the fluctuation-induced in-plane 
conductivity and magnetoconductivity in a series of high quality BaFe2- :r Ni a As2 single 
crystals covering doping levels from the optimal one (x = 0.1) up to the highly overdoped 
(x = 0.2). The sharp superconducting transition allowed us to investigate fluctuation 
effects in almost all the accessible reduced-temperature window above T c . In turn, 
the smooth temperature and field dependences of the normal-state in-plane resistivity 
permitted a reliable estimation of the temperature onset for the fluctuation effects. 
As an example of the usefulness of these data to probe the different approaches for 
the superconducting fluctuations around T c in iron pnictides, we have also presented 
here a detailed comparison with a version of the phenomenological Gaussian Ginzburg- 
Landau approach that includes an energy cutoff, which reduces to the well-known 
Aslamazov-Larking result at low reduced-temperatures and magnetic fields. Our data 
are in good agreement with this approach in all the studied temperature and field 
regions, suggesting the absence, even in the very low magnetic field regime, of local 
superconducting order associated with phase fluctuations, [3] |6| at present a debated 
aspect of the HTSC.[13l [HJ US EH EH] Other contributions to the fluctuation-induced 
in-plane magnetoconductivity (Maki-Thompson, Zeeman, and DOS) are found to be 
negligible, as is also the case in HTSC.[50J Finally, we have also shown that when 
analyzing the low-field regime, it is crucial to take into account the T c -inhomogeneities 
associated with chemical disorder, which are always present to some extent even in the 
best crystals, due to the unavoidable random distribution of doping ions.[2T] It would 
be interesting to extend our present analysis to the critical region by using the LLL 
scaling proposed in Ref. [2], and also check our present results in other families of Fe- 
based superconductors, in particular in the more anisotropic 1111 pnictides, where the 
fluctuations dimensionality is at present a debated issue. [TJ [HJ [9j [TQl [HI [12] 
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contrast with our present results, and in spite of the three-dimensional character of 
this compound, it is claimed that the fluctuation conductivity agrees with the two- 
dimensional Aslamazov-Larkin expression, and also that the onset temperature is as 
high as ~ 2.5T C . 




Figure 1. (a) Detail of the resistive transition corresponding to the x = 0.1 crystal. 
The thick line is the background contribution. Inset: overview up to room temperature 
for all crystals studied, (b) Temperature dependence of dp a b/dT indicating the onset 
of fluctuation effects for the x = 0.1 crystal. The line corresponds to the normal-state 
background, (c) x dependence of T c and T onset . The blue line is a fit of a degree- two 
polynomial to T c (x). The red line, obtained as 1.5T c (x), is in relatively good agreement 
with the observed T onset values. The factor 1.5 is what one would expect from the GGL 
approach presented in Appendix B by using a cutoff constant of c = 0.4. 
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Figure 2. Example (for the x — 0.1 crystal) of the H dependence of the fluctuation- 
induced in-plane magnetoconductivity (a), and conductivity (b), for temperatures just 
above T c . The lines correspond to Aa a b(0) — Aa a b(H) and Aa a b{H) respectively, 
evaluated by using Eq. @ and the same £ a &(0), £ c (0) and c values. The coincidence 
between Aa a b(0) — Aa a b(H) and <7 a f,(0) — cr a b{H) (which is background-independent), 
represents an important verification of the adequacy of the procedure used to subtract 
the background contribution. The dashed line in (b), evaluated by using h = e in 
Eq. @, represents the crossover to the region at which finite field effects are important 
(shadowed). [35] 
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Figure 3. Examples of the Aa a b dependence on e for the crystals with x = 0.1 (a) 
and 0.15 (b). The conventional AL approach [dashed line in (a)] fits the data in the 
low-e and -if region but is not valid above e ~ 0.1. In contrast, Eq. ([2]) (solid lines) 
reproduces the Aa a b reduction at high-e and -if values. The £ab(0) and £ c (0) values 
resulting from the fits of Eq. © for all doping levels are presented in the inset, where 
the error bars include the uncertainty in the p a b amplitude. For comparison, the AL 
approach evaluated with the same parameters is represented as a dotted line in (a). 
The upturns at low- if (indicated by arrows) are due to T c inhomogeneities associated 
with the random distribution of Ni dopants. The dot-dashed line in (b) corresponds 
to the effective-medium approach proposed in Ref. [3D] to account for the effect of T c 
inhomogeneities, evaluated by assuming a Gaussian T c distribution 1.2 K wide. 
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Appendix A. Uncertainty in Aa a b(e, H) associated with the uncertainties in 
the transition temperature and in the normal-state contribution. 

Here we present some representative examples of how the uncertainty in the transition 
temperature and in the normal state contribution affect the determination of the 
fluctuation contribution to the in-plane electrical conductivity. 

In our work T c is estimated as the temperature at which the slope of the p a b{T) 
curve is maximum. In turn, the T c uncertainty was estimated as AT C = 2(T C — T~), 
where T~ is the temperature at which p a b vanishes. In Fig. I Alf a) we present a detail of 
the zero-field resistive transition of the optimally-doped crystal (with x = 0.1) showing 
the locations of T c , T~ , and T c + = T c + AT c /2 (the latter represents an upper limit for 
the T c value). In Fig. IAl( b) we show the Aa ab dependence on the reduced temperature, 
e = ln(T/T c ), as obtained by using those three representative T c values. As may be 
clearly seen, Aa a b is almost unaffected by the uncertainty in T c above e = 0.03. The 
line is the best fit of Eq. ([2]) above this e value to the solid data points. In the presence of 
a finite magnetic field Aa a b presents a smoother behavior at T c (its divergence is shifted 
to lower temperatures) and Aa a b(e) is less affected by the uncertainty in T c . This is 
illustrated in Fig. IAl( c). where the result corresponding to p^H = 5 T is presented. 
The solid line in this figure is the best fit of Eq. ([2]) to the solid data points. 

The same analysis is presented in Fig. IA2I for one of the overdoped samples 
(x = 0.15). In this case, as may be seen in the detail of the resistive transition presented 
in Fig. IA2( a). AT C /T C is larger than in the optimally-doped crystal. As commented 
above, this may be attributed to the random distribution of Ni dopants and to the T c 
dependence on the doping level. This effect is also present in doped cuprates (see, e.g., 
Ref. [21]) and may be intrinsic to non-optimally-doped superconductors. The Aa a b 
dependence on the reduced temperature, as obtained by using the three representative 
T c values shown in Fig. IA2( a). is presented in Fig. IA2( b). Circles and rhombus were 
obtained with p H = 1 T and 5 T, respectively. The effect on Aa a b of changing T c by 
T~ or T c + may still be accounted for by Eq. (T5]) by just changing the £ c (0) value (which 
affects the Ao~ ab amplitude) by ±20%. 

Next, we present an example (corresponding to the x = 0.1 crystal) of the typical 
uncertainty in Aa a b associated with the determination of the background contribution. 
In Fig. IA3( a) we present a detail of the temperature dependence of the in-plane 
resistivity in the normal state up to ~ 4T C . The lines are the background contributions as 
determined by fitting a degree-two polynomial in different temperature intervals between 
the onset temperature for the fluctuation effects (T onset ) and 4T C . The robustness of the 
background contribution to changes in the fitting region is a consequence of the smooth 
behavior of the normal-state in-plane resistivity in a wide temperature region above 
Tonset (if varies about 1% from T c up to ~ 1.5T C ). The same background contributions 
are presented in the dp a b/dT vs. T representation in Fig. IA3( b). The well defined linear 
behavior of dpabjdT up to ~ 4T C justifies the use of a quadratic form to determine 
Pab,B(T). This last figure also illustrates that T onset (estimated as the temperature at 
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which dp a b/dT rises above the extrapolated normal-state behavior beyond the noise 
level) is almost independent of changes in the background fitting region. In Fig. IA3( c) 
we present the ^-dependence of Aa a b as obtained by using the background contributions 
in Fig. IA3( a). Solid lines correspond to Eq. fl2]) evaluated with £ c (0) = 0.8 nm and 
the values indicated for the cutoff constant. As is clearly shown, the uncertainty in 
the background leads to an uncertainty in the cutoff constant below ±10%. Given the 
relation c = ln(T onset /T c ), this leads to an uncertainty in T onset which remains below 
±4%. 
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Figure Al. (a) Detail of the resistive transition corresponding to the x = 0.1 crystal, 
indicating the upper, lower, and midpoint T c values. The effects of the T c choice on 
the resulting Acr a 6(e) is presented in (b) and (c) for hqH = and 5 T, respectively. 
The lines are fits of Eq. ^ to the solid data points. 
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Figure A2. (a) Detail of the resistive transition corresponding to the x = 0.15 crystal, 
indicating the upper, lower, and midpoint T c values. The effects of the T c choice on 
the resulting A<r a b(e) is presented in (b) for fioH — 1 and 5 T. The solid (dashed) lines 
are fits of Eq. ([2]) to the data obtained with T~ (T c + ). The difference in the resulting 
£ c (0) is about ±20%. 
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Figure A3, (a) Detail of the normal-state in-plane resistivity for the x = 0.1 crystal, 
showing the robustness of the background contribution to changes in the temperature 
interval used to determine it. The same backgrounds are presented in (b) in the 
dpab/dT vs. T representation. The effects of the background choice on the resulting 
Aa a f,(e) is presented in (c), where the lines are fits of Eq. @ to the different data sets. 
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Appendix B. 3D anisotropic GGL approach for Aa a b in the finite-field 
regime 

Our starting point is the model proposed by A. Schmid (adapted for a 3D anisotropic 
superconductor), based on a combination of the standard GL-expression for the 
thermally-averaged current density of the superconducting condensate, 



A 2 

' e -(**W- tfW*) — A(t)<**#) 



% m vnrc 
with the generalized Langevin equation of the order parameter, 

irhao 



(B.l) 
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- + — $(r) 

dt he y ' 



a e + 



2m* 



2e A , N 



tt(r,t) + G(r,t). 



(B.2) 



In these equations, m* is the mass of the Cooper pairs, c is the speed of light, $ is the 
electrical potential, k B is the Boltzmann constant, = h 2 /2m*£ 2 (0) is the Ginzburg- 
Landau normalization constant [here £(0) is the GL coherence length amplitude], and A 
is the vector potential. Note also that Eq. flB.2j) almost coincides with the conventional 
time- dependent Ginzburg-Landau equation of the order parameter. The only difference 
is the presence of a random force, G(r,t), which must be completely uncorrelated in 
space and time. The latter implies that G(r,t) will verify 

(G*(r, t)G(r', t')) = a5(r - r')S(t - t'), (B.3) 

where a is a normalization constant that may be determined by just taking into account 
that in the stationary limit Eq. ( IB. 21) has to reproduce the equilibrium thermal average 
of the squared order parameter. |38j This directly leads to a = nhao/i. 

The presence of an homogeneous electrical field, E, applied at the instant t — 
may be taken into account in this formalism by applying 

r-cE ti f t >o 

W \ if t < V ; 

and $(r) =0 to Eqs. fIB.ll) and (IB. 21) . After a standard Fourier-like expansion of the 
order parameter this leads to 



2e 



J = ^ _1 E— (^P + 2eEt)(|M/ p | 2 ), 



(B.5) 



and, respectively, 
dt 
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*p(t) + G p (t), 



(B.6) 
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where ^ p (t) is the Fourier- component of the order parameter corresponding to the 
wavevector p = (p x ,P y ,Pz) and G p (t) represents the random force in momentum space 
that, accordingly to Eq. ( IB. 31) . will verify 

(G* p (t)G p/ (t')} = a5(p - p')5(t - t'). (B.7) 

The thermally averaged current density at an arbitrarily high electrical field may be now 
obtained by introducing in Eq. (IB .511 the (|^/ p | 2 )-expression resulting from Eq. (IB. 61) . 
The latter is a differential equation with solution 

lknT r n r l , 

dtG p (t) x 

(hp + 2eEt"f 



exp 



Trha 
8k B T c0 



Trfiao 



dt" 



a e + 



2m* 



(B.8) 



Then, by using Eqs. (IB.7j) and (IB.8j) . we obtain 
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dt' x 
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(B.9) 



and, subsequently, the supercurrent density at high applied electrical fields will be given 
by 

^(^p + 2eEt) / dt' 



32ek B T^ 
7rhm*aoV ^ 
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exp 
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(hp + 2eEt 
2m* 



ii\2 



(B.10) 



As addressed in Ref. [38], to solve the time- integrations involved in Eq. (IB.lOj) it is 
convenient to introduce first the factor (fop + 2eEt) inside the integral over t' and, then, 
to apply the following changes of variables 

hp = Kk- eE(t + 1') 

t' = u + t (B.ll) 



t" 



u')+t . 



Equation (IB.lOj) is then transformed into 
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an expression that solving the integration over u' and introducing the dimensionless 

reac 

e 2 e(0)nE 



variable x = l&k ? T eu reads as 
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where E* = ~^~|^ e3//2 is a temperature-dependent electrical field characteristic of 
each material that governs the non-Ohmic regime of the fluctuation induced supercurrent 
that appears for E > E* . Note also that in this last expression we have already omitted 
de frk-term that appears in the integral over u in Eq. (IB. 12[) . The reason is that such 
term is odd and, thus, the contributions of the addends corresponding to k and — k 
cancel each other when carrying out the sum over the momentum spectrum. 

In what follows we will restrict ourselves to applied electrical fields verifying 
E «C E* . The last term in the exponential of Eq. (1B.13j) can be then suppressed, and 
the current density exhibits the linear dependence with the electrical field characteristic 
of the Ohmic-regime. By using J = ActE and performing the trivial integral over x, we 
find the following expression for the fluctuation-induced conductivity as sum over the 
modes of the spectrum of the fluctuations 

Aa = ^n^?[i + e( e )k 2 p ' (R14) 

that transforming the /c-sumations into fc-integrals through 

k 

may be rewritten as 

e 2 £ 2 (0) f f f dk x dk y dk z 

J J J [ £ + £2( )k 2 ] 2 ' (R16) 

Equation (lB.16j) has been derived in the framework of a general model for a 3D 
isotropic superconductor. However, in view of the ratio between the superconducting 
coherence length amplitudes in BaFe 2 _ x Ni x .As 2 superconductors, the 3D anisotropic 
scenario seems to be more appropriate for these compounds. In this last dimensional 
case, the scale variation of the order parameter in each spatial direction is determined 
by the corresponding superconducting coherence length. Thus, considering that the x 
and y directions lay on the afe-plane and that z corresponds to the crystallographic 
c-axis, Eq. ( IB. 161) can be adapted to anisotropic 3D superconductors by applying 
£ 2 (0)k 2 ->■ £ 2 fe (0)£; 2 + £ 2 fe (0)£; 2 + £ 2 (0)fc 2 [here £ a &(0) is the in-plane superconducting 
coherence length amplitude]. Using polar coordinates for the xy-plane this leads tcQ 

a e 2 £afe(0) r r k xy dk xy dk z 

ab J J [e+e ab (o)ki y +e c (o)kir [ ■ n 

If an external magnetic field H is applied parallel to the c— direction, the in-plane 
spectrum of the fluctuations becomes equivalent to that of a charged particle in a 
magnetic field. (53] Thus, k xy in Eq. (IB. 171) must be replaced by ie ^ H (n+ |V where fio 

* Note that this equation is different from the one that may be obtained by means of the Kubo formula 
[see, e.g., Ref. [52] . Both integral cores show essentially a fc~ y 3 -dependence, but achieved with different 
fc xy -powers in the numerator and denominator. Thus, at zero applied magnetic field and without limits 
on the fc-integrals both methods lead to the same result [Eq. ((TJ]. However for H ^ or considering 
a cutoff in the energy spectrum the final Ac-expressions show slight quantative (but not qualitative) 
differences that deserve further studies [cf. Ref. [52] and, respectively, Ref. [18].] 
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is the vacuum magnetic permeability and n = 0, 1 . . . is the Landau- level index. As a 
consequence, the integral with respect to k xy is transformed into a sum over n through 
j- J k xy dk xy — > J2n- Besides, we must also include as a multiplier to Eq. ( IB. 17)) the 
so-called Landau degeneracy factor given by = (here O is the magnetic flux 
quantum). The resulting expression for Aa ab is 

A(Jab = Y^hI dkz ^ [e + h[2n + + S(°)*3" a ' (R18) 

where h = H/H c2 (0) is the reduced magnetic field and H c2 (0) = cj) /27rfj, ^ b (0) ^ s the 
upper critical magnetic field perpendicular to the afr-planes, linearly extrapolated to 
T = K. 

To take into account the limits imposed by the uncertainty principle to the shrinkage 
of the superconducting wavefunction when e or h increase, an energy cutoff must be 
applied to Eq. ( 1B.18j) .|[T9] This restricts the sum over n and the integration over k z 
through n max = — 1 and \k™ ax \ = a/c — e/£ c (0) (here c is a cutoff constant expected 
to be of the order of 0.5).[l9~| 134] leading to 
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(B.19) 



In the zero magnetic field limit, i.e., for h e, c, this equation is transformed into 

arctan </— arctan 



Aa. 



ab 



16^ c (0) 



(B.20) 



that corresponds to the paraconductivity under an energy cutoff. At low reduced 
temperatures and magnetic fields, for h,e c, the cutoff effects become 

unimportant. [T9l 131] Accordingly, in this regime Eqs. (IB. 191) and (1B.20I) reduce to the 
c-independent expressions 

Aa ab = — 4 ~, \H [°°dx^ ( £ 4A + x 2 ) (B.21) 



32^< c (0) V h Jo r V 2h 
and, respectively, the AL paraconductivity in a 3D anisotropic superconductor 
[Eq. (fl])].[27] Note finally that Eqs. (IB"T9]) or f lR2tJ]) lead to the Aa ab vanishing at 



c. 
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